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GLUON RESUMMATION IN VECTOR BOSON PRODUCTION AND DECAY 1 

C.-P. YUAN 

Department of Physics and Astronomy, Michigan State University, 
East Lansing, MI 48824, USA 

After a brief review on the technique of resumming the large logarithmic terms a" hi m (Q 2 /Qr) 
due to the effects of multiple gluon emission predicted by the QCD theory, I discuss its 
application in the production of W , Z and 77 in hadron collision. 

1 Introduction 

It is the prediction of the QCD theory that at hadron colliders the production of Drell-Yan pairs 
or weak gauge bosons (W^ and Z) are often accompanied by gluon radiation. Therefore, to 
test the QCD theory or to probe the electroweak properties in the vector boson productions, it 
is necessary to include the effects of multiple gluon emission. At the Tevatron (a pp collider), 
we expect about 2 x 10 6 W ± and 6 x 10 5 Z bosons produced at yfS = 1.8 TeV, per lOOpb -1 
of luminosity. This large sample of data is useful for (i) QCD studies (single and multiple scale 
cases), (ii) precision measurement of the W boson mass and width, and (hi) probing for new 
physics (e.g., Z'), etc. To achieve the above physics goals requires detailed information on the 
distributions of the rapidity and the transverse momentum of W^/Z bosons and of their decay 
products. 

Consider the production process /11/12 — > VX. Denote Qt and Q to be the transverse 
momentum and the invariant mass of the vector boson V, respectively. When Qt ~ Q, there 
is only one hard scale in this problem, and the fixed-order perturbation calculation is reliable. 
When Qt <C Q, this becomes a two scale problem, and the convergence of the conventional 
perturbative expansion becomes impaired. Hence, it is necessary to apply the technique of QCD 
resummation to resum the singular terms: 

^T~^{«s(£ + 1) +a|(L 3 + L 2 ) +a!(L 5 + L 4 )+a 4 5 (L 7 + L 6 ) + ... 

+a 2 (L + 1) +a|(L 3 + L 2 )+a 4 (L 5 + L 4 ) + ... 

+a|(L+ l) + a|(L 3 + L 2 ) + ...}, 

where L denotes ln(Q 2 /Q T ) and the explicit coefficients multiplying the logs are suppressed. 

Resummation of large logarithms yields a Sudakov form factor and cures divergence at 
Qt — * 0. It was pioneered by Dokshitzer, D'yakonov and Troyan (DDT) who performed the 
analysis in Qy-space which lead to the leading-log resummation formalism a. Later, Parisi- 
Petronzio showed! that for large Q the Qt — ► region can be calculated perturbatively by 
imposing the condition of the transverse momentum conservation, 
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in the fe-space (the impact parameter, the Fourier conjugate to Qt)- This improved formalism 
sums also some sub leading- logs. As Q — > oo, events at Qt ~ may be obtained asymptotically 
by the emission of at least two gluons whose transverse momenta are not small and add to zero. 
The intercept at Qt = is predicted to bei 



da 



dQl 



A 2 
Q 2 



no 



where r/ = Aln 1 + \ with A = 12C F /(33 - 2n f ), and r] ~ 0.6 for n f = 4 and C F = 4/3. 

Collins-Soper-Sterman (CSS) extended I the work by Parisi-Petronzio in 6-space and applied 
the properties of the renormalization group invariance to set up a formalism that resums all the 
large log terms to all orders in a s . This is the formalism I will concentrate on in this talk. 

Recently, there is another theory model in Q^-space (extension of DDT) proposed by Ellis 
and Veseliu, which does not have either the exact transverse momentum conservation or the 
renormalization group invariance conditions. The reader can find a detailed discussion of this 
formalism in the talk by K. Ellis at this conference. Despite of the imperfectness of the formalism, 
it can still be useful for the and Z physics program at the Tevatron. I shall come back to 
this point in the conclusion section. 



2 Collins-Soper-Sterman (CSS) resummation formalism 

In the resummation formalism by Collins, Soper and Stermani, the cross section is written 
in the form 



da(h x h 2 -> VX) 1 



dQ 2 dQ\dy {lir f 



d 2 b e^W(b, Q,x l ,x 2 ) + Y(Q T , Q, xi , x 2 ) 



where Y is the regular piece which can be obtained by subtracting the singular terms from the 
exact fixed-order result. W satisfies a renormalization group equation. Its solution is of the 
form 

W(Q, b, Xl ,x 2 ) = e" S(gACllC # (Jr^, b, xx,x^j , 
where the Sudakov exponent is defined as 



f c iQ 2 du 2 
S(Q,b,C 1 ,C 2 )= -4- 



A (a s (ji), d) In ( j +B (a s (p), C ls C 2 ) 



and the x\ and x 2 dependence of W factorizes as 

W {§tb ,b,Xl,X2 ) = ^>^ {§tb ,b,Xl ) Cjh2 {^tb AX2 ) ' 



where Cjh is a convolution of the parton distribution with a calculable Wilson coefficient, called 
Cj a function: 



Cjh(Q, b,x)=^2 j — C ja 



fx C 3 ( 

0,fl= -TiQ) Ja/h C,A« 



X 



b 



where a sums over incoming partons, and j denotes the quark flavors with (electroweak) charge 
ej. A few comments about this formalism is listed below: 
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• The A, B and C functions can be calculated order-by-order in a s . 

• A special choice of the renormalization constants can be made so that the singular terms 
obtained from the expansion of the CSS resummed calculation agrees with that from the 
fixed-order calculation. This is the canonical choice. It has C\ = C3 = 2e~ 7B = 60 and 
C2 = C1/60 = 1, where 7^ = 0.577. . . is Euler's constant. 

• b is integrated from to 00. For b S> 1/Aqcd, the perturbative calculation is no longer 
reliable. Hence, a non-perturbative function is needed in this formalism to compare theory 
prediction with experimental data. 

We refer the readers to Ref. i for a more detailed discussion on how to apply the CSS resum- 
mation formalism to phenomenological physics at hadron colliders. 



2. 1 Non-perturbative function 

As noted in the previous section, it is necessary to include a non-perturbative function in the 
CSS resummation formalism to incorporate some longdistance physics that is not accounted by 
the perturbative derivation. Collins-Soper postulated^ 



W fk (b) = W fk (K)Wf k P (b), 



with 



VI + Jb/br. 



so that b never exceeds b max and W^ip*) can be reliably calculated perturbatively. (In a 
numerical calculation, b max is set to be, say, 1/(2 GeV).) Based upon a renormalization group 
analysis, they found that the non-perturbative function can be generally written as 



Wjjfc (b, Q,Qo,xi,x 2 ) = exp 



Fj/hifa,b) - F k/h 2 ( x 2,b) 



where Fx, Fj/ hl and F^, h2 must be extracted from data with the constraint 



Wjj- p (b = 0) 



1. 



Furthermore, F\ only depends on Q, while Fj/ hl and F^, h2 in general depend on x\ or x 2 - Later, 
in Ref. i, it was shown that the F\(p)\n(%x) dependence is also suggested by analyzing the 



infrared renormalon contribution in Qt distribution. 



2. 2 Testing the universality of W^ p 

The CSS resummation formalism suggests the non-perturbative function to be universal. Its role 
is similar to the parton distribution function (PDF) in any fixed order perturbative calculation, 
and its value needs to be determined by data. The first attempt to determine such a universal 
non-perturbative function was done by Davies, Webber and Stirling (DWS)Ejn 1985 using Duke 
and Owens parton distribution function. In 1994, Ladinsky and Yuan (LY) El observed that the 
prediction of the DWS set of W% p largely deviates from the R209 data (p + p — > + X 

at y/S = 62 GeV) using the CTEQ2M PDF. To incorporate possible ln(r) dependence, LY 
postulated 



W jjc (b,Q,Q ,xi,x 2 ) = exp 



-9ib 2 



g 2 b 2 In 



_Q_ 
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Figure 1: Comparison of 4pb CDF-Z data with two different theory model predictions. 



where x x x 2 = r. A "two-stage fit" of the R209, CDF-Z (4pb~* data) and E288 (p + Cu) data 
gave gi = 0.11±g;^ GeV 2 , g 2 = 0.58j^ GeV 2 , and g 3 = -1.5lg;J GeV" 1 , for Q = 1.6 GeV and 
frmaz = 0.5 GeV -1 . Unfortunately, a fortran code error in calculating the parton densities inside 
the neutron lead to a wrong g$ value. (An erratum with corrected values will be submitted 
after the completion of our revised analysis.) Currently, Brock, Ladinsky, Landry, and Yuan El 
are revisiting this problem using the R209 (p + p), CDF-Z (p + p with 4pb _1 data), E288 
(p + Cu), and E605 (p + Cu) data, with CTEQ3M PDF. The preliminary results show that, 
for Qo = 1.6 GeV and b max = 1/(2 GeV), two forms for W% p (b, Q, Qo, x\, x 2 ) give good fits 
(x 2 /dof ~ 1.4): (i) g\ = 0.24, g 2 = 0.34 and 53 = 0.0 (DWS form, pure Gaussian form in 
6-space, without x dependence), and (ii) g\ = 0.15, g 2 = 0.48 and 53 = —0.58 (LY form, 
with a linear b term and x dependence). We are in the process of determining the uncertainties 
of these fitted parameters 

2.3 Run- IB W/Z data at the Tevatron 

The Run-IB W/Z data at the Tevatron can be useful to test the universality and the x depen- 
dence of the non-perturbative function W^ p (b, Q, Qo, x±, x 2 ). In Fig. p], we show the prediction 
of the two different global fits (2-parameter and 3-parameter fits) obtained in the previous section 
using the CTEQ3M PDF. (The CTEQ4M PDF gives the similar results.) 

We note that for Qt > 10 GeV, the non-perturbative function has little effect on the Qt 
distribution although in principle it affects the whole Qt range (up to Q). This is also clearly 
shown in the figures of Ref. i. (Because of the limited space of this article, I will not reproduce 
here the figures that can be found in Ref. i.) To study the resolution power of the Tevatron 
Run-IB Z data on determining the non-perturbative function, we have performed a "toy global 
fit" as follows. First, we generate a set of "fake Run-IB Z data" (assuming 5,000 reconstructed 
Z bosons) using the original LY fit (g\ = 0.11, g 2 = 0.58 and 53 = —1.5). Then, we combine this 
set of "fake" data with the low energy Drell-Yan data as listed above to perform a global fit. 
Using the 3-parameter form (the LY form) , we get back the g\ and g 2 values from the fit but the 
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gs value is smaller by a factor of 2. The best fit gives g\ = 0.11, <?2 = 0.57 and gz = —0.88. (This 
amounts to shift the prediction on the mass and the width of the W boson by about 5 MeV and 
lOMeV, respectively 113.) With this large sample of Z data, it can be clearly illustrated that a 
single parameter without Q dependence (i.e. the parameter g± alone) cannot yield a good global 
fit. 

3 W/Z Production and Decay at the 1.8 TeV Tevatron 

In Ref. i, we have presented the results of a detailed study on the distributions of the decay 
leptons from the W/Z boson produced at the Tevatron. Here, we shall only present a few of the 
most interesting results. 

3.1 "Matching" from low to high Qt 

First of all, let us examine the distribution of the transverse momentum of W boson produced at 
the 1.8 TeV Tevatron. It is obvious that the CSS formalism gives automatic matching in the Qt 
distribution if A, B, C and Y functions were known to all orders in a s . However, in practice we 
only calculate those functions to some fixed order in a s , therefore, some "switching" procedure 
(to switch from the total resummed result to the fixed-order perturbative result) should be 
applied to have a continuous distribution that would match the fixed-order result in the large 
transverse momentum Qt (of the order Q) where ln(Q 2 /Q T ) is small and the resummed result 
is less accurate. Our matching prescription is to switch from the resummed prediction to the 
fixed-order perturbative calculation as they cross around Qt ~ Q- This switching procedure is 
done for any given Q (invariant mass) and y (rapidity) of the W/Z boson. We note that this 
procedure is different from that proposed in Ref. E-l As shown in the figures of Ref. B, including 
a higher order Y term, as expected, improves the "matching" between the resummed and the 
fixed-order results. At 0(a s ), that crossing occurs at about 50GeV, while at 0(a 2 ), it occurs 
at about 70 GeV for W boson production. 

3.2 Total cross section 

It was shown by an analytical calculation in Ref.i that if "matching" (switching) is chosen to be 
at Qt = Q, then the CSS resummed rate is the same as the NLO rate. However, in that case, the 
Qt distribution is not smoothly continuous for Qt close to Q. Since the region of Qt > 50 GeV 
only contributes to the total rate by less than 2%, and the CSS resummed formalism includes 
also some even higher order contributions beyond NLO, we regard the difference between the 
resummed rate and the NLO rate as a gauge of the uncalculated higher order contribution which 
is shown to be small at the Tevatron. 



Total cross sections of pp — » (H^ + or Z)X at the Tevatron, calculated in different prescriptions, 

in units of nb. 



V 


E 

J -'cin 


Fixed Order 


Res. (1,1,1) 


Res. (2,1,1) 


Res. (2,1,2) 


Experiment 




(TeV) 


0(a° s ) 0(a s ) 


Pert. 0(as) 


© Pert. 0{as) 


© Pert. 0(a%) 




w+ 


1.8 


8.81 11.1 


11.3 


11.3 


11.4 


11.5 ± 0.7 


w+ 


2.0 


9.71 12.5 


12.6 


12.6 


12.7 




z° 


1.8 


5.23 6.69 


6.79 


6.79 


6.82 


6.86 ± 0.36 


z° 


2.0 


6.11 7.47 


7.52 


7.52 


7.57 





In the above Table, the values of the strong coupling constants used with the CTEQ4L and 
CTEQ4M PDF's are a { g ] {M z ) = 0.132 and a { g ] {M z ) = 0.116, respectively. (Res. (2,1,2) 
denotes the result with A and B calculated to a 2 s order, C to as, and Y to a 2 s order, etc.) 
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pr (GeV) 

Figure 2: NLO (dashed) and resummed G(as) (solid) calculations Using the LY fit with cuts (Qt < 

30 GeV,p T + ' v > 25 GeV). 

We note that the total rate is not sensitive to either non-perturbative function or matching 
prescription. 

3.3 Lepton transverse momentum 

At the Run-2 of Tevatron, the number of interactions per crossing increases, and the transverse 
mass measurement becomes less accurate. It is useful to also measure the inclusive lepton P T 
spectrum. In Fig. ^, we show that the difference between the resummed and the NLO results is 
much larger than the dependence on the non-perturbative function in the resummed calculation 
(although its dependence is not negligible). 

3.4 Lepton charge asymmetry 

The distribution of the rapidity of the charged lepton from W boson decay provides useful 
information to determine the ratio of the up and down quark parton distribution functions. It 
was shown in Ref. i that without imposing the kinematic cuts, the prediction of the resummed 
calculation is the same as that of the NLO. This is obvious because an integration over the 
entire Qt space in the CSS formalism reproduces the NLO rate as discussed above. However, 
with kinematic cuts, the rapidity distribution of the charged lepton or the W boson is differ- 
ent between the resummed and the NLO calculations. Furthermore, for this distribution, the 
dependence on the non-perturbative function is negligible. 

The resummed calculation also predicts different correlations between the two decay leptons 
from Z boson. A couple of examples are the distribution showing the balance in transverse 
momentum Apx = \PtLt \Pt\i an< ^ the angular correlation z = —p T x • p T 2 /[max(p^ , p^)] 2 . We 
refer the reader to Ref. B for more details. 
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yy Resummation at Vs=1.8 TeV 




Figure 3: The predicted distribution for the invariant mass of the photon pair M 77 from the resummed calculation 
compared to the CDF data, with the CDF cuts imposed in the calculation. 



4 Di-photon rates and distributions 

It is straightforward to extend the above CSS formalism to calculate the distribution of photons 
from hih,2 — ► 77X. In Ref. 113, we have included the full content of NLO contributions from qq 
and qg subprocesses. The accuracy of the resummed result is similar to that for the Drell-Yan 
and W/Z productions. In addition, we also included part of the higher order contribution from 
the gg process whose lowest order contribution comes from one-loop box diagrams. This was 
done by including A^2) term in the Sudakov factor to resum higher order large logs due to initial 
state radiation, and the NLO contribution, of 0(a^ m a^), to the hard scattering subprocess is 
added in an approximate fashion. 

We found that the gg — * 77X rate is not small at the Tevatron. Including part of the NLO 
99 ~* 77 approximately doubles the LO gg box contribution to the cross section. This is clearly 
seen in the distribution of the invariant mass of the photon pair, as shown in Fig. ||. The 
other distributions, such as the distribution of the transverse momentum of the photon pair Qt, 
compared with CDF and DO data, can be found in Ref. c3. We note that to accurately predict 
the production rates of the di-photon pairs at the CERN Large Hadron Collider (LHC) requires 
a better calculation of the gg fusion rate. Namely, a full NLO result is needed to predict the 
distribution of the di-photon pair with large Qt or small A0. 

In Ref. 0, we also presented our results for the di-photon production rate and distributions 
at the fixed-target experiment (E706 at Fermilab): p + Be — ► 77 + X at v / S'=31.5 GeV. In this 
case, we found that the gg rate is small, and the distribution of the transverse momentum of the 
photon pair Qt is dominated by the non-perturbative contribution. Therefore, this data can be 
useful for determining the non-perturbative function associated with the gg initiated processes 
in the CSS formalism. We refer the readers to Ref. El for more details. 



5 Conclusion 

In conclusion, the effects of QCD gluon resummation are important in many precision mea- 
surements^ A Monte Carlo package ResBos ( ( Res ummed Bos on Production and Decay)) is 
available t3 for studying the effects of gluon resummation (with NLO hard part corrections) in 
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hadron collisions. 

Before closing, I would like to comment on the Qt space resummation formalism proposed 
in Ref. □, and compare that with the b space resummation formalism (CSS formalism) discussed 
in this article. Despite of its imperfectness in the theory structure, in practice, the Qt space 
resummation formalism may still prove to be useful because it does not require a "switching" in 
the large Qt region and its calculation takes much less CPU time due to the fact that there is 
no need to perform a Fourier transformation from the b space to the Qt space. This is obvious 
for Qt above 10 GeV where the non-perturbative contribution is not important, as discussed 
above. However, in the small Qt region, it is not obvious that both formalisms will give the same 
prediction of the Qt distribution of the boson after the non-perturbative part is fit by the Qt 
of the Z boson using the Tevatron Run-IB and Run-2 data. Nevertheless, if one is not interested 
in testing the universality of the non-perturbative function, as suggested by the CSS formalism, 
then it seems likely that it is possible to choose a proper form of the non-perturbative function 
in the Qt space resummation formalism to reproduce (within the experimental uncertainties) 
the prediction of the b space formalism when considering the Tevatron W/Z data alone. If one 
is interested in the universality property of the formalism]] so that one can use the same values 
of the non-perturbative function to predict future data at different hadron colliders (such as the 
LHC), then it becomes unclear whether these two formalisms will always give the same physics 
prediction. 
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